Enhanced Fermi surface nesting in superconducting BaFe 2 (Asi_. r P a ;)2 revealed by de 

Haas-van Alphen effect 
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The three-dimensional Fermi surface morphology of superconducting BaFe2(Aso.37Po.63)2 with 
Tc — 9K, is determined using the de Haas- van Alphen effect (dHvA). The inner electron pocket 
has a similar area and k z interplane warping to the observed hole pocket, revealing that the Fermi 
surfaces are geometrically well nested in the (rr, 7r) direction. These results are in stark contrast to 
the Fermiology of the non-superconducting phosphides (x — 1) , and therefore suggests an important 
role for nesting in pnictide superconductivity. 



A ubiquitous ingredient in theories of high- 
temperature superconductivity is the role of Fermi 
surface (FS) nesting - a property originating from the 
shape of the FS which enhances quasiparticle scattering 
along particular directions in momentum space. FSs 
that support nesting are typically low dimensional with 
a weakly varying morphology in at least one direction 
in momentum space. In the case of the pnictides, it is 
nesting in the (tt, it) direction between electron and hole 
FSs which is most often invoked [lHEJ . From the earliest 
band structure calculations on the '1111' compounds 
e.g. LaFeAsO, it was clear that this was a possible 
mechanism to explain both the magnetism and the 
superconductivity [5[. However, the '122' pnictides e.g. 
BaFe2As2 do not have such a strongly nested FS, yet 
order equally strongly in the absence of doping 3] . Fur- 
thermore LaFePO, has a FS which is measured to have 
well matched electron and hole pockets, and is therefore 
close to satisfying a nesting condition, yet has no known 
magnetic ordering |6j, [7|. This has led several authors to 
invoke the local magnetic environment, particularly of 
the As-Fe-As tetrahedron as a route to understanding 
these cases @, 01 Q- The role of nesting, is therefore 
highly debated. In order to determine whether a given 
FS supports nesting, an intimate knowledge of the full 
three-dimensional morphology is required and this is the 
purpose of the present study. 

We investigate single crystals of superconducting 
BaFe2(Aso.37Po.63)2, T c — 9K. The compounds 
BaFe2As2 and BaFe2P2 are isoelectronic. The former 
orders antiferromagnetically at ~ 138K, whereas the lat- 
ter is an ordinary metal down to the lowest tempera- 
tures measured. For intermediate compositions there is 
a superconducting dome that peaks at x ~ 0.3, with 
T c = 30K. These materials are particularly useful for 
dHvA studies because the disorder is weak and long mean 
free paths can be maintained. As a consequence, quan- 
tum oscillatory phenomena can be measured, as done 
recently by Shishido et al. Q . This study made two im- 



portant observations: first, the size of the Fermi surface 
increases with x beyond optimal As/P substitution and 
second, the effective mass decreases. However, due to 
the resolution limits of the data the details of the three- 
dimensional curvature of the Fermi surface could not be 
determined for the superconducting compounds and no 
hole pockets were observed. In the present paper we re- 
solve the morphology of the electron pockets and that of 
one hole pocket. We conclude that a large amount of the 
FS is geometrically nested. The details revealed in this 
study provide compelling evidence for the role of nesting 
as a mechanism for high temperature superconductivity 
in this family of compounds. 

High quality single crystals of BaFc2(Asi_ :E P ;E )2 with 
residual in-plane resistivity ratios p(300K)/p(1.8K) ^20, 
were grown from an FeAs flux. The composition x = 0.63 
was determined by microprobe analysis. Torque mag- 
netometry was performed using piezoresistive microcan- 
tilevers in fields of up to 35T in DC Bitter magnets at the 
NHMFL, Florida. The oscillatory part of the torque sig- 
nal T oc M x B originates from the dHvA effect. Each 
dHvA frequency is related to an extremal cross sectional 
area Ak of the FS in momentum space via the Onsager 
relation ¥—{h/2 , ne)Ak. For simply warped, quasi-two di- 
mensional cylinders two frequencies are expected, a maxi- 
mum ('belly' orbit) and a minimum ('neck' orbit). As the 
angle of the field with respect to the sample is changed, 
these extremal orbits traverse larger sections of the Fermi 
surface and at particular angles, known as Yamaji angles, 
all cross-sectional orbits match, leading to an enhanced 
amplitude in the oscillation [10]. Band structure calcula- 
tions were performed for the end members BaFe2As2 and 
BaFe2P2 using an augmented plane wave plus local or- 



bital method as implemented in the WIEN2K code [11 
The composition x = 0.63 does not have a structural 
transition at low temperature and therefore the (room 
temperature) experimentally determined tetragonal unit 
cell parameters are used. These are a = 3.96(3.84) A, 
c = 13.039(12.44) A and z P = 0.3538(0.3456) for the 
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FIG. 1: The background subtracted raw data (a-c) for three different angles taken at T — 0.3K, and corresponding Fourier 
spectra in (d-f). The Fourier amplitudes of (e) and (f) are normalized with respect to (d) for clarity. The 7 pockets have an 
enhance amplitude near 9 « 51° probably due to a Yamaji effect, (g) gives an example of raw torque data and (h) gives the 
Fourier amplitude as a function of temperature with LK fits for the f3 (blue), a (green) and 7 (red) pockets, from which we 
determine the effective mass, (g) and (h) are taken at 9 — 16°. 



As(P) variants. 

Fig. QJg) shows the the raw torque signal for a single 
field sweep up to 35 T. dHvA oscillations are observed 
for fields above ~ 20 T. In Fig. Q] (a-c) we show the 
background subtracted data at three typical angles. The 
Fourier spectrum is generated by implementing discrete 
Fourier transform methods on the signal (in inverse field) 
for each angle and is shown in adjacent panels. At most 
angles the spectrum is dominated by two (split) peaks 
which we label ax, 2 and /3x,2- In total five frequencies are 
observed, which we denote as a%, a2, pi, $2 and 7. We 
observe a further low frequency labeled e which has some 
scatter about a value of 250T at most angles. However, 
its value is sensitive to the background subtraction and 
the field window used and so it is likely a noise artifact 
of the Fourier transform. 

The calculated FS for BaFe2P2 and BaFe2As2 (in the 
tetragonal structure) is shown in Figure [5] (d) and (e) re- 
spectively. Both consist of two similarly warped concen- 
tric electron cylinders at the M point. BaFe2As2 consists 
of three hole sheets while BaFe2P2 consists of only two 
hole pockets centered at the T point. The hole sheets of 
BaFe2P2 are much more warped than BaFe2As2, consis- 
tent with dHvA measurements on SrFe2P2[l3i suggest- 
ing the latter is closer to satisfying a nesting condition. 
The dHvA angle dependence associated with these FSs 
is shown in the Figure [2][a) and (c) for BaFe2P2 and 
BaFe2As2 respectively. 

Figure [2] (b) shows the angle dependence of each of the 
observed frequencies, broadly corresponding to the calcu- 
lated band structure of the undoped compounds, though 
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FIG. 2: (color online), (a) The dHvA angle dependence ex- 
pected from LDA for BaFe2P2- (b) the measured angle depen- 
dence for a superconducting sample of BaFe2(Aso.37Po.63)2- 
(c) The dHvA angle dependence expected from LDA for 
BaFe2As2. The angle 9 = refers to field parallel to the [001] 
crystallographic direction, and 9 = 90 to the [100] (9 = 90) 
crystallographic direction, (d) and (e) show the FS as cal- 
culated for BaFe2P2 and BaFe2As2, respectively. The hole 
pockets are centered at the F point in the zone and the elec- 
tron pockets at the M points. 
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substantial adjustment is required (see below). It should 
be noted that this discrepancy between the experimen- 
tally determined Fermi surface cross-sections and those 
calculated for the end members (x = 0&1) can not be 
reconciled by using the experimentally determined lat- 
tice parameters for the x = 0.63 compound, which only 
weakly affect the calculation. The strongest signal comes 
from the frequencies labeled a and j5. In both LaFePO 0] 
and SrFe2P2 [l2| the strongest signal was from the elec- 
trons pockets, which tend to have the longest mean free 
paths. It is natural to assume that this remains the case 
in the doped compounds and we assign these frequencies 
to the electron pockets. 

The morphology of the pockets is best revealed by plot- 
ting the dHvA frequencies multiplied by cos#, as shown 
in Figure [3] (a). For a two dimensional FS, the frequency 
dependence appears flat, whereas for warped quasi- two- 
dimensional cylinders the extremal cross-sections weave 
between the Yamaji angles. The /3 pocket appears signifi- 
cantly warped (a difference in neck and belly frequencies 
of 175T). a is observable up to 60°, suggesting a sin- 
gle pocket exists around 950T. The assignment oli/ol-2 is 
given to the neck/belly of this pocket and this is con- 
firmed by the comparable effective masses (see Table IJ). 
This suggests that the pocket has a weak warping corre- 
sponding to a neck/belly difference of only 85T. The high 
angle data is vital for the assignment of the 7 pocket. 
Approaching 50° the a frequency abruptly turns down 
(see arrow in Figure [3] (a)), only to appear again near 
its average value a still higher angles. Simultaneously, 
the 7 frequency turns upward, joining the a downturn at 
around 51°. This suggests that for much of the angular 
range the a spectrum in fact contains a third frequency, 
the belly of the 7 pocket. Further confirmation of this 
assignment is evident in the 51° spectrum where the am- 
plitude of the 7 is enhanced beyond that of the a fre- 
quency, indicating the first Yamaji angle for this pocket 
(see Figure [1] (f)). While it is possible that the 7 orbit 
originates from the same Fermi surface as the a orbit, 
this seems unlikely given firstly, that the a orbit remains 
at the average value of ~ 950T at angles > 51° and the 
effective mass of 7 is around twice that of a (see below). 
We conclude that 7 is a hole pocket. Even without fur- 
ther analysis, it is clear that this compound has at least 
one hole pocket which is geometrically well matched to 
the inner electron pocket. 

As all frequencies deviate substantially from calcula- 
tions, we perform rigid band shifts on both the BaFe2As2 
and BaFe2P2 Wien2K calculated band structure, as sum- 
marized in Figure[3] (b) and (c). The (3 pocket can be 
reproduced by the outer electron pocket of either band 
structure with relatively moderate shifts in the Fermi en- 
ergy (see caption Figure [3|) . The average frequency of 
the a pocket can also be reproduced by the inner elec- 
tron pocket of either, but the morphology seems quite 
different from observations in both calculations. The 
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FIG. 3: (color online) (a) The observed dHvA frequencies 
plotted as FcosO. The observed downturn in the a pocket 
joining the Yamaji angle of the 7 pocket is shown, (b) The 
best fitting rigid band shifts to the BaFe2As2 band structure 
with Wien2K orbit numbers labeled. The /3 pocket can be 
broadly matched with a shift of -51meV of band 5, the a 
pocket can be matched with a shift of -71meV of band 4 and 
finally, the 7 pocket can be well approximated with a shift of 
+25meV to band 2, even giving a similar warping to that of 
the data, (c) Similar shifts to the band structure of BaFe2?2 
can give good approximation of the measured /3 pocket, with 
a -106meV shift to band 4, a with a shift of -130meV to 
band 3. However, though a shift of +84meV to band 1 can 
bring the average frequency into alignment with 7, the pocket 
is significantly more warped. Larger shifts cause band 1 to 
form a closed 3D FS. In (a) we also show a fit of a simple 
cosine warped cylinder to the 7 pocket and good agreement 
is achieved with the data, particularly the Yamaji angle. 



band structure in BaFe2As2 and BaFe2P2 predicts sig- 
nificantly more warping than is observed, (see Figure [3] 
(b) and (c)). Similarly fitting the 7 pocket, only the 
the BaFe2As2 band structure fits the data adequately, 
but again the morphology is off because the Yamaji an- 
gle deviates significantly from the data (by about 15°). 
We finally fit the 7 frequency to simple (cosine) warped 
cylinder. Much better agreement is achieved with the 
experimental data and the Yamaji angle at ~ 51° is al- 
most perfectly matched. In summary, the experimentally 
determined FS can support nesting in a manner that is 
not obvious from band structure calculations of the end 
members. 

The dHvA mean free path and effective mass are ex- 
tracted in the conventional manner using the Lifshitz- 
Kosevich (LK) formula [l3| and the results are summa- 
rized in Table |U As expected, the hole pocket has a sig- 
nificantly shorter mean free path than the electron pock- 
ets. For the effective mass, the a and /? pockets are in 
good agreement with those reported in RefH The effec- 
tive mass of 7 has the largest error bar because of the 
spectral leakage from low frequency noise in the Fourier 
spectrum, which tends to enhance the mass. Neverthe- 
less, it is clear that the effective mass is around twice that 
of the electron pockets. In order to determine the renor- 
malization we compare our masses to that of the shifted 
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FIG. 4: (color online). Two dimensional sections through the 
FS in the (110) plane shown for unshifted band structures 
shown in (a) for BaFe2As2 and (c) for BaFe2P2- In (b) we 
show the electron pockets after shifting the BaFe2As2 bands 
to agree with our dHvA, and the a cosine warped cylinder is 
shown for the hole pocket. 



band structure of BaFc2As2 in Table HI We could have 
equally chosen to compare to BaFe2P2, in which case the 
renormalization would be roughly double that shown in 
the table. The mass enhancements are different for elec- 
tron and hole pockets, being about A = m* /mi, — 1 ~ 0.8 
for the a and (3 pockets and A = m* /nib — 1 ~ 2.6 for the 
7 orbit. These enhancements are much larger than ex- 
pected from electron-phonon coupling alone(A ep ~ 0.25) 
[14j . Interestingly, the non-superconducting compounds 
have electron pockets with a greater enhancement than 

This study 



the holes, in contrast to the present case [12 1 
illustrates that as the superconducting dome is entered 
this electron/hole mass asymmetry is reversed, suggest- 
ing a vital role for many body effects. 

The total electron count from a and f3 is 0.1 elec- 
trons per unit cell, whereas the the 7 pocket accounts 
for only .03 holes per unit cell (approximately equiva- 
lent to the a electron pocket), leaving .07 uncompensated 
holes. As such one further hole pocket must exist that 
is not presently observed. Intriguingly, the contribution 
of the a and 7 pocket almost perfectly compensate each 
other, and so the remaining FS must compensate the /3 
pocket. This suggests the intriguing possibility that the 
remaining FS is of a similar size to the f3 pocket, which 
may therefore sustain FS nesting between the two outer 
electron/hole pockets, in addition to the inner. 

dHvA studies of the non-superconduct ing '122' ar- 
senides (x = 0) reveal a reconstructed FS[15j, whereas 
in the phosphides (x = 1) they reveal an unrescon- 
structed FS but which is far from satisfying a nesting 
condition [l2| . In stark contrast, the present work shows 
that the intermediate superconducting compositions have 
well-matched electron and hole FSs, and are therefore 
close to nesting. The differences between our experimen- 
tal determination and the band structure calculation for 
each end member are summarized in Figure H Further- 
more, the superconducting compounds have significantly 
smaller FSs than expected by band structure, in line with 
inter-band scattering theory of Ortenzi et al. between 
well-nested Fermi surfaces [lfj. At smaller values of x 
with a higher T c it is likely that this trend continues and 



TABLE I: Measured dHvA frequencies, effective masses (m*) 
and mean free paths (I), along with the values from the band 
structure calculations. The experimental masses were deter- 
mined at 9 = 16°. The band structure values are also es- 
timated at the same angle from band structure calculations. 
Mean-free-paths could not be calculated separately for both 
neck and belly orbits due to the data resolution. 

Experiment Calculations 
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the hole pockets even more closely matched the electron 
pockets, becoming ever more difficult to observe in dHvA 
studies and perhaps explaining their absence in the data 
of RefE 

In conclusion, we have measured the FS of super- 
conducting BaFe2(Aso.37Po.63)2 mapping the morphol- 
ogy and mass enhancement in detail for both electron 
pockets and one hole pocket. We find that the warping 
of the inner electron FS is more two-dimensional than 
expected in the unshifted band structure, and shares a 
similar degree of warping to that of the hole pocket. The 
correlation between favorable nesting and superconduc- 
tivity presently observed in BaFe2(Asi_ a; P a ;)2 provides 
strong evidence for the importance of nesting in under- 
standing iron-pnictide superconductivity, as suggested by 
a number of theoretical treatments [l|-|5[. 
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